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Abstract 

We investigate two cosmological axion problems ( isocurvature fluctuations and 
domain-wall formation ) in chaotic inflationary universe. It is believed that these 
problems are solved if potential for the Peccei-Quinn scalar field is very flat. How- 
ever, we find that too many domain walls are produced through parametric reso- 
nance decay of the Peccei-Quinn scalar field. Only the axion model with N = 1(N: 
QCD anomaly factor) is consistent with observations. We also point out that the 
flat potential is naturally obtained in a supersymmetric extension of the Peccei- 
Quinn model. If Peccei-Quinn breaking scale F a is about 10 12 GeV, this model 
predicts anisotropies of cosmic microwave background radiation due to the axion 
isocurvature fluctuations which may be detectable in future observations. 
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1 Introduction 



The axion |]T], [3], is a Nambu-Goldstone boson associated with breaking of the Peccei- 
Quinn symmetry which was invented as a solution to the strong CP problem in QCD ||. 
The Peccei-Quinn symmetry breaking scale F a is stringently constrained by laboratory 
experiments, astrophysics and cosmology: allowed range of F a lies between 10 10 GeV and 
10 12 GeV ||. The axion is also cosmologically attractive since it can be a cold dark matter 
if F a takes the higher values ~ 10 12 GeV. 

Inflationary universe || was proposed to solve various problems in the standard 
cosmology (e.g. flatness and horizon problems). There have been constructed many 
inflation models since then. The chaotic inflation model || seems the most attractive 
candidate among them that realizes naturally the inflationary universe. However, when we 
consider the axion in the chaotic inflationary universe, we are confronted with two serious 
cosmological problems, both of which are associated with large quantum fluctuations 
generated in the exponentially expanding phase (i.e., inflationary epoch) of the early 
universe. One is the domain wall problem [|10J. In the inflationary universe fluctuations 
of the axion field a(x) are given by 5a = H/(2tt) where H is the Hubble constant. Since 
a phase 9 a of the Peccei-Quinn scalar field is related to the axion field a(x) by 9 a = a/F a , 
the fluctuations of 9 a is given by 

m 



2nF a 

In the chaotic inflation model, H ~ 10 14 GeV |l l] is required to explain the anisotropics 
of cosmic microwave background radiation (CBR) observed by COBE Then, the 



fluctuations of the phase 9 a become 0(1) for F a < 10 12 GeV, which means that the phase 
is quite random during the inflation. Therefore, when the universe cools down to about 
1 GeV and the axion potential is formed, the axion sits at different positions of the 
potential in different regions of the universe. Since the axion potential has N discrete 
minima (N: QCD anomaly factor), domain walls are produced pi. The domain wall with 



iV = 1 is a disk-like object whose boundary is an axionic string and it collapses quickly 
due to its surface tension. Thus the domain wall with iV = 1 is cosmologically harmless. 
However, the domain wall with N > 2 is disastrous because it forms a complicated network 
with an axionic string and dominates the energy density of the universe quickly [|14[ . 



The second problem is that the quantum fluctuations for the axion field cause too 



large anisotropies of CBR [15, 16, 17, 18, 19]. Since the axion does not have a potential 



during the inflation, the axion fluctuations 5a do not contribute to the energy density 
of the universe. In this sense, the axion fluctuations are isocurvature. After the axion 
acquires a mass m a , its fluctuations become density fluctuations given by 5p a /p a ~ 59 a /9 a 
which causes the CBR temperature fluctuations 5T/T ~ 59 a /9 a . From eq.([TJ), we see that 
the produced CBR anisotropies are 0(1) which contradicts the observation. 



It has been pointed out in ref. JT7, THJ that these problems are simultaneously solved 
if the potential for the Peccei-Quinn scalar field is very flat. In this case it is natural that 
the Peccei-Quinn scalar field $ takes a very large value ~ M p \ during the inflation. Here, 
M p i is the Planck mass (M p \ = 1.2 x 10 19 GeV). The axion field a(x) is defined as 

$(x) = 0(x)exp(m(x)/|(0)|), (2) 
4>{x) : real, — oo < (j) < oo, 
a(x) : real, - 7r/2 < a(x)/\(<f>)\ < 7r/2. 

In the true vacuum <fr takes \(<f))\ = F a . But the (f) takes most likely \{<p)\ ~ M p \ at the 
inflationary epoch. From eq.(|l]) it is clear that the phase fluctuations 59 a are suppressed 
as 59 ' a ~ H/(2irM p \). Therefore the isocurvature fluctuations are suppressed by a factor 
F a /M p \. In ref. |T9| , cosmological effects of the axionic isocurvature density fluctuations 
have been investigated in details by calculating spectra for matter fluctuations and 5T/T, 
and it has been found that the effects of the isocurvature fluctuations are consistent with 
observations in a large parameter region. 

As for the domain wall problem, it is naively expected that the same mechanism 



that suppresses the isocurvature fluctuations does also suppress the production of domain 
walls: i.e. the large effective "F " = | (0) (inflation ~ M p \ during the inflation almost fixes 
the phase of $ which reduces the do main- wall production rate to negligible amount. In 
fact, the $ begins to oscillate after the inflation with its phase being fixed. However, 
the homogeneous Peccei-Quinn field decays into axions as well as <fi particles through 
parametric resonances. The emitted axions induce large fluctuations of the axion field 
(i.e. the phase 56 a ) and produce a large number of domain walls. 

In this letter we study the above cosmological axion problems ( isocurvature-fluctuation 
and domain-wall problems ) in chaotic inflationary universe. We show that the domain 
walls are still produced through the parametric resonance decay even if we take a very 
flat potential for the Peccei-Quinn scalar. The domain wall problem is only avoided for 
the model with N — 1. Furthermore, we also show that the flat potential for $ required 
to solve the isocurvature-fluctuation problem is naturally obtained in a supersymmetric 
(SUSY) extension of the Peccei-Quinn model. Namely, the presence of flat directions is 
a generic feature in SUSY theories and the masses of the fields corresponding to the flat 
directions only comes from soft SUSY breaking terms which are of order the weak scale. 
Therefore the flat potential for $ is naturally realized in the framework of SUSY. 

2 Domain Wall Problem for a Quartic Potential 

Let us first consider the case for a scalar potential given by 



where we have assumed that the potential for the inflaton x is ^X 4 /4- Isocurvature 
fluctuations for cosmologically relevant scales are generated when \ — 4M p i. Therefore 




(3) 



with g <C 1. In the inflationary epoch, the slowly evolves as 




(4) 
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the CBR anisotropies from the axion isocurvature fluctuations is given by [] 

t Pa V37r y \M pl e a J e a 1 ; 

Using a COBE constraint 5T/T < 1(T 5 , we get g< lO -11 for # a ~ 0(1). Notice that the 
coupling constant A for the inflaton is already very small as A ~ 10~ 13 |Ll| . 

Next we consider domain wall production with the potential ([3]). The Peccei-Quinn 
field $ starts to oscillate after the inflation. We only consider the evolution of $ in 
the "radial" direction (j){x) with its phase being fixed since the quantum fluctuations 
of the axion phase 89 a = 5a/ (0) are negligiblly small when |(0)| ~ M p \ as discussed 
in the introduction. As easily seen from eq.(|3]) the "radial" potential has two minima 
(0 = ±F a ). Since the initial potential energy density of is large, the oscillates 
beyond the potential hill at = 0. The loses its energy through cosmic expansion and 
particle production and settles down to one of the potential minima eventually. 

If the dynamics after inflation is purely classical, the final value of (F a or —F a ) 
depends on the initial value of 0. However, it has been recently pointed out that the 
oscillating coherent field quickly losses its energy through violent particle production 
due to parametric resonance P0 |. The parametric resonance decay occurs during the 



oscillation. In the present model, only a self-production (e.g. (0) + (0) — » $ + $) is 
effectivefj (Notice that the "radial" oscillation induces excitations of the phase (=axions) 
as well as 0.) In this case the classical picture does not apply at all. The important 

fact is that the quantum fluctuations due to the emitted and axion particles becomes 

^or a more accurate analysis including comparison with the most recent observations, see ref. 
2 Even if one assumes the inflaton potential ~ m^x 2 /2, one gets essentially the same conclusion as in 

the text. 

3 Note that in the conventional notation of the axion field $ = K^le^W/IWI the two minima corre- 
spond to 6® and 9® + ir where 9° a is the initial value of the axion phase. 

4 ln principle, the $ can couple to other fields ?y. Although the mass of 77 is large ~ 0(F a ) in the 

true vacuum of the coherent $ field also decays into 77 particles through the parametric resonance. 
However, this does not change the conclusion in this paper as will be noted in the end of this section. 



quite large. [] Contributions of the quantum fluctuations 5(f) and 5a = \<p\56 a to the 
energy density of $ are ~ g(5cj)) 2 M^ l |2(| and ~ g(69 a ) 2 M^ h f\ respectively Since most 
of the initial oscillation energy gM^ is transferred to $ particles through the parametric 
resonance, we expect that the fluctuations of <fi and 9 a are given by 5<fi ~ M p \ and 59 a ~ 1, 
respectively, which leads to production of domain walls. 

For more quantitative analysis, we need to solve equations for modes of the fluctuations 
5<pk and 5a k which correspond to <p particle and axion with momentum k, respectively. 
However it is more convenient to use two field X and Y which are defined by X = i?e($) 
and Y = im($) with initial condition X(0) = |0(O)| and Y(0) = 0. Then the fluctuations 
of <p and 9 a are approximately given by 5<f) ~ 5X and 56 a ~ 5Y/\X\. The evolutions of 
X, 5X and 5Y are described by 

X + gX{X 2 -F 2 a ) = -3g((5X) 2 )X-g((5Y) 2 )X, (6) 
5X k + [k 2 -gF 2 + 3gX 2 ]5X k = -3g((5X) 2 )5X k - g((5Y) 2 )5X k , (7) 
5Y k + [k 2 -gF 2 + gX 2 ]5Y k = -g((5X) 2 )5Y k - 3g((5Y) 2 )5Y k , (8) 

where we have used the mean field approximation ((5X) 3 ~ 35X ((5X) 2 ) ■ •) and ne- 
glected the cosmological expansion.^] Before we present results of the full numerical inte- 
gration of the above equations, we briefly discuss instability of the equations for 5X and 
5Y. When the amplitude of the X oscillation is much larger than F a , a solution to eq.(Q) 
is approximately given by 

X ~ Zsm(c^Zt), (9) 
5 Notice that the wavelengths of the fluctuations induced by the decay are much shorter than the hori- 
zon scale. Therefore, these fluctuations do not make any contribution to the fluctuations with wavelengths 

relevant to large scale structures of the universe. 

6 This energy density comes from the kinetic term of the axion field. 

7 The cosmological expansion reduces the effect of parametric resonance in general. However we find 
that the reduction effect is not large enough to avoid the conclusion in this paper. 
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where c = 0(1) and Z is the amplitude of the oscillation. Using this solution and 
neglecting back reactions, eqs.(|7|) and (§) become well-known Mathieu equation: 

5X'^ + [A x (k)-2q x cos(2z)}SX k ~ 0, (10) 
5Y k " + [A Y {k) — 2q Y cos(2z)]SY k ~ 0, (11) 

where A x = ^"ff + 2<?x, qx = A Y = ^z0- + 2gy, q Y = 5?, z = c^JgZt and 
a prime denotes a derivative with respect to z. For gx,y ~ 1, the Mathieu equation 
has instability for k which satisfies A(k) ~ 1,4, 9, Since A x (k) > 1 the strongest 
instability (resonance) for SX k occurs in the second instability band, i.e. A x (k) ~ 4. On 
the other hand the resonance for SYk occurs in the first instability band and the instability 
is stronger than that for 8X k - Since SYj. roughly corresponds to the axion with momentum 
k it is expected that the homogeneous X field decays dominantly into axions. However, 
as the fluctuations increase, back reactions, i.e., terms in RHS of eqs.(^) and (|8|), become 
significant and the instability is weakened. To see the evolution of the fluctuations with 
back reactions, we need numerical calculations as shown below. 

Precise time evolution of the fluctuations is obtained by full numerical integration of 
eqs. (H) - (H) and results are shown in Fig. [I] - [§ As seen in Fig. |] and ||, the fluctuations 
of both X and Y increase through the parametric resonances. Since the instability of 
5Y is stronger, the fluctuations of Y increase faster than those of X. As the fluctuations 
increase, the back reactions become important and weaken the instability and finally 
the fluctuations become almost constant. As is expected, the final fluctuations are large 
and, in particular, the phase fluctuations are of order 1 (59 a ~ 1), which results in the 
production of too many domain walls (see Fig. |j). 

The above mechanism for producing domain walls seems quite general. In fact we 
have solved equations for the fluctuations with a potential V = g/M^Q^l 2 — F%) 3 and 
obtain almost the same result. We also consider interaction between $ and other scalar 
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fields rj: d n t = — f\&\ 2 \v\ 2 -^ K the coupling / is larger than g, the $ field first decays into 
^-particles. However, the back reactions of 77-particles suppress the resonance decay when 
5r] ~ M p i |2(J . Then the remaining $ field may decay into 0-particles and axions through 
the parametric resonance, which results in large phase fluctuations (~ 0(1)). Thus the 
production of 77 particles cannot suppress the violent fluctuations of the axion phase. 
Therefore, the production of domain walls may not be avoided in the chaotic inflationary 
universe. Since the produced domain walls dominate the density of the universe for the 
axion model with N > 2, only the model with N = 1 is allowed. Fortunately, the N — 1 
model is easily constructed in a hadronic axion model for example pi . 



3 Supersymmetric Potential 

We now consider a SUSY extension of a hadronic axion model. We take the following 
simple superpotential: 

W = h^ + V_-F 2 )V , (12) 

where and \l/o are chiral superfields with Peccei-Quinn charges +1,-1 and 0, 

respectively, and the coupling constant h is assumed to be ~ 0(1).[] Then the scalar 
potential is written as 

Vsusy = /> 2 |$+$- - F 2 a \ 2 + h 2 (\$ + \ 2 + |$-| 2 )|$or (13) 

where $ + , $_ and $0 are scalar components of superfields and ^0, respectively. 
Here we should note that there are flat directions that satisfy 

= F 2 a} $ = 0. (14) 

8 Possible trilinear coupling contributing to $ decay is Ci„t — —f'Qi] 2 . This interaction, however, 

gives unbounded potential at the vacuum |$| = F a and hence we discard this possibility. 

9 We also assume a pair of quark Q and antiquark Q whose Peccei-Quinn charges are —1/2. This 

charges are chosen so that Q and Q have a Yukawa coupling W — JQQ^>+. This model has the QCD 

anomaly factor N = 1. 
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We also add soft SUSY breaking terms to the potential. Then, the scalar potential V 
that we should study is 

V = Vsvsy + <|$ + | 2 + m 2 _|$_| 2 + m 2 |$ | 2 , (15) 

where m± t o are soft masses of O(lOOGeV). 

We discuss cosmological evolution of these scalar fields. The chaotic situation at the 
Planck time may lead to large initial values of order M p \ for $ + , $_ and $o- The soft 
masses can be negligible at the inflation epoch since they are much smaller than the 
expansion rate of the universe (H ~ 10 14 GeV). Therefore these fields roll down into 
the flat direction given by eq . (|14|) . However, one scalar field, either $ + or may stay 
at the initial position, i.e., |$+| or |$_| ~ M v \. We assume that it is the case and take, 
for definiteness, |<&+| ~ M p \ and eliminate $_ and $o using eq.(|TJ|). Then, the effective 
lagrangian for $ + is written as 

L = L kin -V, (16) 

V = m^+l 2 + m 2 F 4 |$+|- 2 . (18) 

Since the potential ( p~8|) is very steep at $ + = and the kinetic term is always positive, 
the energy conservation never allows the $ + to pass through the origin ($ + = 0). Thus, 
this potential has effectively only one minimum. From the above lagrangian we obtain 
equations of motion for the homogeneous field $+ and phase fluctuations SO: 



10 At the inflation epoch, non-vanishing vacuum energy density breaks SUSY and may generally give 
scalar masses of the order of the Hubble constant H p2[ . However, there exist a class of supergravity 
models where such scalar masses are negligible compared with H |23| . We assume such a class of models 
in this paper. 
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where the prime denotes d/(m+dt), (p = \<& + \/F a and 9% ~ / dx 3 58e tkx . In eg. (|I~9|) the 
back reactions of 56 are assumed to be neglected. Although eq.(^) looks different from 
the Mathieu equation, numerical calculations show that it has indeed strong instability 
(see Fig. [5]). Therefore we expect that the parametric resonance induces large phase 
fluctuations and results in the production of domain walls. However these domain walls 
are cosmologically harmless since we consider the hadronic axion model with N = 1 as 
discussed in the previous section. 

Let us now estimate isocurvature fluctuations of the axion in the SUSY model. First 
we rewrite <3>+ and $_ as 

Q-|- 

$-1- = t> ± expi^= , — oo < v± < oo. (21) 

V2\v±\ 



We also define the fields a and b by 



V r«-, (22) 



( v 2 + + y 2_y/2 + {v 2 + + vlY 



V. 



+ 

" ~ (vl+V*)V* a++ (Vl + ^)l/2 fl - ( 23 ) 

From eq.(|13|), the potential V& for b is given by 

V b = -2h 2 F 2 a v + v^ cos + v -) l ' 2 h ^ (24) 

Thus the mass of b is ~ (y\ + v 2 _ ) 1 ^ 2 ~ M v \ ~ 10 19 GeV, which is much larger than the 
Hubble constant during the inflation and hence b quickly becomes zero at the inflation 
epoch. On the other hand, the potential for a is flat and is regarded as the axion. At the 
inflation epoch, the quantum fluctuations of a are given by, assuming v + 3> V-, 

5a ~ 5a + ~ — . (25) 
2ir 

Since the mass of b is very large compared with H, the quantum fluctuations of b are 
negligible and hence <5a_ ~ — (v-/v + )8a+. Therefore, the fluctuations for phases (9±) of 
$± are given by 

S9 ± = ^ - . (26) 

9 



Since v + ~ M p \ and H ~ 10 GeV, fluctuations of the phase 8± are quite small (58± ~ 
10~ 6 — 10~ 5 ). After the axion obtains a mass, the phase fluctuations result in isocurvature 
density fluctuations and produce the CBR anisotropics as explained in the non-SUSY case. 
When $ + and settle down to the true minimum of the potential (i.e., v + ~ t>_ ~ F a ), 
the axion field is written as a ~ a + — a_ and its fluctuation 5a is given by 5a ~ F a 56 + , 
which gives 5p a /p a ~ 59 + /9 + . Thus the 5T/T is predicted as 10~ 6 — 10~ 5 , which is 
consistent with the COBE observation. In particular, if F a ~ 10 12 GeV, the predicted 
CBR anisotropies due to the isocurvature density fluctuations are large enough to be 
detected by future observations as stressed in ref. (19| . 

In SUSY Peccei-Quinn models, there exists an axino which is a supersymmetric partner 
of the axion. The mass of the axino is presumed to be the SUSY scale, i.e. 100 GeV 
- 1 TeV. The axino of such a large mass is cosmologically harmless since it decays into 
lighter particles quickly. 



4 Conclusion 

In summary, we have investigated two cosmological axion problems (isocurvature fluc- 
tuations and domain-wall formation) in chaotic inflationary universe. The isocurvature- 
fluctuation problem is solved if potential for the Peccei-Quinn scalar field is very flat. 
However, we have found that too many domain walls are produced through the paramet- 
ric resonance decay of the Peccei-Quinn scalar field and only the model with N = 1 is 
allowed. We have also pointed out that the flat potential which is necessary to solve the 
isocurvature-fluctuation problem is naturally obtained in a SUSY extension of the Peccei- 
Quinn model. Namely, we do not need any small coupling constant in the SUSY model 
for obtaining such a flat potential. On the other hand a very small coupling constant 
(g ^ 10~ n ) is necessary in non SUSY Peccei-Quinn models. 
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Figure Captions 



Fig.l Time evolution of the fluctuations of X. z is the time in units of (c v /^Z(0)) _1 
where c is the 0(1) constant, g is the self-coupling constant and Z(0) is the initial 
amplitude of the oscillation of the Peccei-Quinn scalar field. We take X(0) = M p \/3, 
and g = 1CT 13 . 

Fig. 2 Time evolution of the fluctuations of 5X. 

Fig. 3 Time evolution of the fluctuations of SY. 

Fig. 4 Time evolution of the fluctuations of 9 a .We define 59 as arctan(y / (5Y) 2 /\X\). 

Fig. 5 Time evolution of the phase fluctuations with k = 0.5m + . We take |$ + (0)| = 
lOF a and m = m + = m_. This figure shows clearly the instability of the phase 
fluctuations. 
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Figure 1: Time evolution of the fluctuations of X. z is the time in units of (c v /^Z(0))~ 1 
where c is the 0(1) constant, g is the self-coupling constant and Z(0) is the initial am- 
plitude of the oscillation of the Peccei-Quinn scalar field. We take X(0) = M p i/3, and 
g = 1(T 13 . 
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Figure 2: Time evolution of the fluctuations of 8X. 
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Figure 5: Time evolution of the phase fluctuations with k = 0.5m + . We take |<3> + (0)| = 
10F a and m = m + = ra_. This figure shows clearly the instability of the phase fluctua- 
tions. 
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